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Accurate estimation of third-order moments from turbulence 
measurements 
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Abstract. Politano and Pouquet's law, a generalization 
of Kolmogorov's four-fifths law to incompressible MHD, 
makes it possible to measure the energy cascade rate in in- 
compressible MHD turbulence by means of third-order mo- 
ments. In hydrodynamics, accurate measurement of third- 
order moments requires large amounts of data because the 
probability distributions of velocity-differences are nearly 
symmetric and the third-order moments are relatively small. 
Measurements of the energy cascade rate in solar wind tur- 
bulence have recently been performed for the first time, but 
without careful consideration of the accuracy or statistical 
uncertainty of the required third-order moments. This pa- 
per investigates the statistical convergence of third-order mo- 
ments as a function of the sample size N. It is shown that 
the accuracy of the third-moment ((<5v||) 3 ) depends on the 
number of correlation lengths spanned by the data set and a 
method of estimating the statistical uncertainty of the third- 
moment is developed. The technique is illustrated using both 
wind tunnel data and solar wind data. 



1 Introduction 

In the solar wind, coupling between large- and small-scale 
turbulence occurs at kinetic scales defined by the ion gyro- 
radius and the ion gyro-period. At these scales, the turbulent 
energy cascade undergoes a transition from large magneto- 
hydrodynamic (MHD) scales to small plasma kinetic scales 
where the energy is ultimately dissipated by collisionless 
processes. Detailed understanding of the energy cascade pro- 
cess at MHD-scales is a prerequisite for studies of this cou- 
pling. Here we focus on one particular aspect of MHD-scale 
turbulence which is of some practical importance, namely, 
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the determination of the energy cascade rate from measured 
data. 

MHD-scale turbulence in the solar wind is often mod- 
eled using the theory of incompressible MHD because of 
its relative simplicity, even though the solar wind is known 
to be compressible. In the solar wind, the energy density 
of MHD turbule nce is comparable t o the plasma thermal 
energy at 1 AU dBelcher and Davisl 1 1 97 lb and the turbu- 
lent energy cascade is believed to significantly heat the so- 
lar wind plasma as it flows from ~ 1 AU to several tens of 
AU. Theoretical work has shown that plasma heating caused 
by dissipation of the turbulence can likely explain the ob- 
served radial temperature profile of the solar wind which de- 
creases more slo wly than would be the c a se if the expansion 



were adiabatic dMatthaeus et all 1996; ZanketaU 11999; 



ic (Mattha 
all Il999t 



Matth aeus et all Il999t ISmith et all I200U llsenberg et al 
2003b . To refine these theories, accurate measurements of 



the energy cascade rate are needed. Recently, the energy 
cascade rate e has been directly measured for the first time 
in the solar wi nd using a gene r alization of Kolmogorov's 
four-fifths law (jMacBride et all 120051 ISorriso-Valvo et al. . 
120071 iMacBride et all 120081; iMarino et all 12008b . Before 
discussing this, it may be helpful to provide some back- 
ground information on Kolmogorov's four-fifths law. 

For turbulent flows in ordinary incompressible fluids such 
as air or water the energy cascade rate e is often measured 
indirectly by means of the energy dissipation rate 
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dv 

dx 



(1) 



where v is the kinematic viscosity and the coefficient 15 
arises from the assumption that the turbulence is isotropic 
dPopel l2000l p. 134). The energy cascade rate can also be 
measured directly by means of Kolmogorov's four-fifths law 
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valid for isotropic turbulence, where 

Svu (r) = [v(x + r) — v(x)] ■ e r 



(3) 



is the component of the velocity fluctuation in the direction 
of the disp lacement r and the lengfhscale r lies in the iner- 
tial range dKolmogorovi 1 1 99 ll: iFrischi 1 1 995b . Note that Kol- 
mogorov's four-fifths law (fJI is independent of the kinematic 
viscosity v and can be applied even when the kinematic vis- 
cosity is unknown, but the accurate evaluation of the third- 
order moment (12) requires much more data than the second- 
order moment ([TJ. 

Kolmogorov's four-fifths law was originally derived 
for homogeneous isotropic turbulence and a similar 
law was later derived by Monin for hom ogeneous 
anisotro pic turbulence; see IPodesta et al.l (2007) for refer- 
ences. Politano and Pouquet ( 1998a b) generalized these 
fundamental results of Kolmogorov and Monin from the the- 
ory of incompressible hydrodynamic turbulence to incom- 
pressible MHD turbulence. It is important to emphasize 
that Politano and Pouquet's law holds for both isotropic and 
anisotropic turbu lence, although this fact was not explic- 
itly mentioned by Politano and Pouquet ( 1998ah . This is es- 
pecially important in MHD where statistical isotropy may 
not hold in the presence of an ambient magnetic field. A 
derivation of Politano and Pouquet's law which is similar to 
Fr isch's de r ivati on of Kolmogorov's four-fifths law is given 
bv |Podestal(l2008h . 

Politano and Pouquet's law has recently been ap- 
plied to obtain direct measurements of the energy cas- 
cade rate in the solar wind under the simplifying as 
sump ti on that the turbulen c e is isotropic (j MacBrid e et al 
20051: ISorriso-Valyo et all l2007t [MacBride et all 12008 



Marino et al.L 120081) . iMacBride et all (12008b have also inves- 
tigated a non-isotropic 1D/2D hybrid model that is believed 
to be descriptive of the solar wind. The method used in all 
these studies consists of the evaluation of certain third-order 
moments which are similar to those in equation (J2J, except 
that for incompressible MHD turbulence the relevant third- 
order moments contain combinations of velocity and mag- 
netic field fluctuations (or, equivalently, fluctuations in the 
Elsasser variables). From the linear scaling of these third- 
order moments, the energy cascade rate is obtained without 
any knowledge of the dissipation processes or the viscous 
and resistive dissipation coefficients in the the solar wind. 

The solar wind studies mentioned above have not given 
careful consideration to the convergence properties of third- 
order moments which raises the question: how much data 
is required to accurately estimate the third-order moments? 
The study by Sorriso-Valvo et al. (2007) used approximately 
2000 data points to compute the third-order moments while 
the study by MacBride et al. (2008) used close to 10 6 data 
points. The purpose of the present work is to investigate the 
accuracy of third-order moments as a function of the sample 



size N (the number of data points used in the analysis). An 
important conclusion is that the accuracy of third-order mo- 
ments depends on the number of correlation lengths spanned 
by the data set (defined below). The number of correlation 
lengths determines the accuracy and statistical uncertainty of 
third-order moments computed from measured data, not the 
number of data points N. It turns out that for turbulence stud- 
ies where the skewness of the distribution is usually small 
the accurate estimation of third-order moments requires large 
amounts of data. The reason is partly because the third mo- 
ment is not an absolute moment (|a;| 3 ) but a signed moment 
(a; 3 ) and, therefore, is subject to cancellation effects. The 
theory describing the convergenge of these moments is illus- 
trated using turbulence data from the ONERA/Modane wind 
tunnel. The same techniques can be applied to third-order 
moments in solar wind turbulence which exhibit similar be- 
havior. 



2 Theory 

2.1 Uncorrelated time series 



Given N independent samples xi,X2,.. 
domly from a probability distribution f(x) 
the distribution can be estimated as 
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etc. Now focus attention on the third moment and let 

N 



(5) 



(6) 



(7) 



n=l 



Note that Als(N) is itself a random variable whose probabil- 
ity distribution can be derived, in principle, from the pdf f(x) 
of the random variable x. Now suppose that we know the 
mean and standard deviation of the random variable M^(N) 
denoted by ^3 and 03, respectively. If (a; 3 ) 7^ 0, then for 
M-s(N) to be an accurate estimate of (x 3 ) the standard devi- 
ation must be small compared to the mean, that is, 



a 3 (N) 



< 1. 



(8) 



This condition can be used to estimate the sample size N 
required to obtain an accurate estimate of the third moment. 
Hereafter, it is assumed that (a; 3 ) 7^ 0. 
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From equation (0 and the independence of the x n , the first 
and second moments of A/3 (N) are 



1 N 



(9) 



and 



N N 
n— 1 m=l 

1 , «, 



N 



TV- 1 
TV 
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(10) 



Thus, the variance is 



a 2 (N) = < [M 3 (TV)] 2 ) <A/ 3 (TV)) 2 = - [(x 6 ) 



and 



M3 



r 3\2 



- 1 



1/2 



(ID 

(12) 



The value of TV required to make the last equation small (<C 
1) depends on the ratio (x 6 )/(x 3 ) 2 and, therefore, depends 
on the distribution function f(x). If the ratio (x 6 )/(x 3 ) 2 is 
on the order of unity, then TV > 10 3 may be adequate. But, 
if this ratio is much larger than unity, then TV will have to in- 
crease accordingly. The relation (fT2l) shows that to increase 
the accuracy of the third-moment by a factor of ten requires 
an increase in the sample size TV by a factor of 100. This 
slow rate of convergence imposes practical limitations on es- 
timates of third-order moments obtained from experimental 
data. 



2.2 Correlated time series 

For applications to turbulence, the random variable x 
is a velocity difference such as <5u|| and the sequence 
Xi,%2, ■ ■ ■ , xn is usually not mutually stochastically inde- 
pendent. For example, two velocity increments that over- 
lap in space or time are usually correlated to some degree. 
In this case, the number of "independent samples" TV in the 
above theory should be replaced by the number of correlation 
lengths of the quantity under consideration. For the third- 
order moment A/3 it is necessary to use the correlation length 
or correlation time r c of the time series y n = x^ and replace 
TV in equation (12) by the number of correlation lengths 



T 

N c = -, 



(13) 



where T — Nt s is the temporal record length and t s is the 
sampling time. Thus, equation (12) takes the modified form 



03 

/'3 



1 



N, 



(x 3 ) 2 



1/2 



(14) 



or, equivalently, 



1/2 



r3\2 



- 1 



1/2 



(15) 



where r c = nt s . Note that this is almost the same as equation 
(12) except for an additional scale factor n 1 / 2 . Because the 



correlation times of the sequences x n and y n 



can be 



different it is important to use the correlation time t c of the 
sequence y n = x 3 in equations (fT3ll-(fT5]l when analyzing 
the third-order moment. 



3 Textbook example 

An example is now given to illustrate the theory described 
in section 2. Consider the slightly skewed distribution func- 
tion 

f (x) = -^=[l + a{x + a)]e-^ 2 '\ (16) 



2tt 

where the parameter a characterizes the skewness of the dis- 
tribution. The distribution has zero mean, (x) = 0, and re- 
duces to the Gaussian distribution when a = 0. If y = x + a, 
the pdf of y is 



f(y) 



1 



'2tt 



(1 + ay)e- y2/2 . 



(17) 



The moments of the distribution function ( fT71 ) can be com- 
puted from the characteristic function 



F(k) 



f(y)e iky dy = (l + iak)e- k ^ 2 (18) 



by means of the well known relations 

F'(0)=i(y), 
F"(0) = t 2 (y 2 ), 
F"'(0) = i 3 (y 3 ), 



(19) 



etc. After some tedious calculations, the first six moments 
are found to be 

(y) = a, 
(y 2 ) = 1, 



<2/ 3 > 

(y 4 ) 
(y 5 ) 
(y 6 ) 

Thus, from the relation (x n ) 
ments of the variable x are 

(x) = 0, 
(x 2 ) = 1-a 2 , 
(x 3 ) = 2a 3 , 
(x 4 



3a, 
3, 

15a, 
15. 



(20) 



((y — a)"), the first six mo- 



3a , 



(21) 



(x 5 ) = 20a 3 + 4a 5 , 
(x 6 ) = 15 -45a 2 -45a 4 



5a b 
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For the particular distribution function ( fT6] l, the ratio ( fT2] > 
takes the form 



CT3 


1 


15 - 45a 2 - 45a 4 - 5a 6 


M3 




4a 6 



1/2 



For a = 0.25, for example, this becomes 



1± 

M3 



111 

Vn' 



(22) 



(23) 



This ratio is small if N > 10 6 . In this idealized example 
where the distribution function f(x) is known, the number of 
samples required to obtain an accurate estimate of the third 
moment from experimental data can be computed explicitly. 
For turbulence data acquired in the laboratory, such precise 
estimates cannot be computed a priori because the distribu- 
tion function f(x) is unknown. 



4 A practical approach 

When working with experimental turbulence data the dis- 
tribution function f(x) is usually unknown so that the ratio 
(x 6 )/(x 3 ) 2 in equation ( fl~2l cannot be evaluated. A prac- 
tical approach is to compute the third-moment M3 (N) from 
the data and then construct the empirical distribution function 
for A/3 (A), where N is now fixed (a constant). This too may 
be impractical because of the large number of data points re- 
quired. However, the distribution function for M3 (A) con- 
tains more information than is needed. Just a few indepen- 
dent estimates of the third-moment M3 (A), perhaps 10, may 
be sufficient to obtain a rough estimate of the ratio in equa- 
tion (0. The number of samples A can then be increased 
until the ratio so obtained satisfies the inequality ([SJ. This 
is a simple way of controlling the accuracy of third-order 
moments estimated from turbulence measurements. The ef- 
fectiveness of the method can be improved by increasing the 
number of independent estimates of M3 (A) used to compute 
the mean and standard deviation. The standard deviation ob- 
tained from the data provides a rough estimate of the 1-cr 
error for the third-order moment. 

A more precise analysis can be performed by comput- 
ing histograms, means, and standard deviations of the third- 
moment M3 (A) for progressively larger values of A. Accu- 
rate values of the mean /i3(A) and standard deviation 03(A) 
can be obtained for values of A much smaller than the record 
size. Fitting the measured ratio I03//Z3I to the functional 
form A/^fN, where A is an adjustable parameter, it is then 
possible to extrapolate the ratio | cr 3 / /x 3 1 to larger A where 
direct calculations have poor statistics or are unattainable as 
a consequence of the limited record size. An alternate fit- 
ting function is A/N p where A and p are two adjustable 
parameters. From this extrapolation it is possible to deter- 
mine the value of the sample size A required to achieve any 
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Fig. 1. Power spectrum of longitudinal velocity fluctuations mea- 
sured in the Modane wind tunnel (blue dots). For comparison, the 
red line is proportional to / _6 , Kolmogorov's theoretical inertial 
range spectrum. The inertial range extends from approximately 10 
Hz to 10 3 Hz. The onset of the dissipation range is indicated by the 
change in slope around 10 3 Hz. 

desired accuracy of the ratio | cr 3 / /j, 3 1 and, therefore, of the 
third-order moment A/3 (A). This approach is accurate and 
effective as long as sufficient data are available and requires 
no apriori knowledge of the distribution function f(x) or its 
moments. The same technique can also be applied to accu- 
rately determine moments of any order provided sufficient 
data are available. 



5 Illustration using wind tunnel data 

The technique described in the previous section shall 
now be applied to study turbulence data from the ONERA 
wind tunnel in Modane, France, ch aracterized by a Taylor- 
scale Reynolds number R\ 
Malecot et al 



2500(KahalerrasetaL, 1998 



2000; Ga gne et aUl2004l) . This particular data 
set is a time series consisting of 1.44 x 10 7 data points with 
a sampling rate of 25 kHz and an average velocity of 20.37 
m/s. The inertial range extends from ~ 10 Hz to ~ 10 3 Hz as 
inferred from the power spectrum shown in Figure 1 . Now, 
consider the third-order moment 



((^||) 3 >=(Ki)-^ + T)] 3 ) 



(24) 



where r = 1/50 s or, equivalently, / = 1/r = 50 Hz. This 
time lag is chosen for study because it lies inside the inertial 
range displayed in Figure 1 . 
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The third-order moment is computed using a contiguous 
series of N data points. A set of N contiguous data points is 
called a data block. A series of successive data blocks are 
then used to compute a series of third-order moments, one 
for each data block. The first data point in a given data block 
is separated from the first point of the next successive data 
block by an offset m where, ideally, m — N. When the sam- 
ple size N is not small compared to the record size, smaller 
values of the offset m are used so that the total number of data 
blocks yields a sufficient statistical sample. Note, however, 
that when the offset m is smaller than N the third-order mo- 
ments obtained from successive data blocks become depen- 
dent (because the blocks overlap) and, consequently, to ob- 
tain good statistics it is advisable not to let m become much 
smaller than N. This tradeoff is unavoidable when working 
with records of finite length. 

The set of third-order moments obtained for a given sam- 
ple size N are used to generate a histogram of third-order 
moments as shown in Figure 2. The number of third-order 
moments N mom is equal to the number of data blocks em- 
ployed in the calculation. As expected, the width of the dis- 
tributions as measured by the standard deviation is a decreas- 
ing function of N. Moreover, the results for the ratio of the 
standard deviation to the mean are in approximate agreement 
with the N~ x / 2 convergence rate predicted by the theory in 
section 2. Thus, increasing N by a factor of ten causes a de- 
crease in the ratio of the standard deviation to the mean by a 
factor of ~ 3. 

The third-order moment obtained using all available data 
is M 3 = —0.406 (m/s) 2 . Note that the mean /J3 displayed 
in Figure 2 is approximately independent of N. This is to be 
expected because for any sequence of numbers partioned into 
successive non-overlapping blocks, the average of the mean 
values for each data block is equal to the mean value of the 
entire record. In practice, the data blocks may not completely 
cover the given record because the record size is not divisible 
by N and, therefore, the equality is only approximate. This 
explains the approximate independence of n 3 {N) versus N 
in Figure 2. See also equations (O and ( fTTl i which predict 
that ti 3 (N) is independent of N and a 3 (N) ~ 1/N 1 / 2 . 

How much data is required to obtain an accurate estimate 
of the third-order moment? This depends, of course, on 
the level of error which is tolerable for the application at 
hand. The relative error is measured by the ratio |<t 3 / 1 . 
This quantity is plotted as a function of N in Figure 3 (up- 
per plot). To ensure a reasonably large number of third- 
moments N mom , the offset m between successive data blocks 
IS 772 = N when N < 10 5 and 777 = N/2 when N > 10 5 . 
The number of third-moments N mom is shown in the middle 
plot in Figure 3. The theoretical relation (15) takes the form 



M3(iV) 



304 

y/N' 



(25) 



where the value 304 is obtained using the empirical values 
of the sixth moment Mq = 169.5 (m/s) 6 , the third moment 
M 3 = -0.4062 (m/s) 6 , and the correlation time r c ~ 90t; s 
defined as the time where the autocorrelation function equals 
1/2 (Figure 3). Inspection of the theoretical curve, the red 
line in Figure 3, shows that to achieve the level of precision 
I03/V3I < 0.1 would require N > 10 7 data points or, equiv- 
alentTy/V c > l.lx 10 5 correlation lengths. The total number 
of data points contained in the data set is 1.44 x 10 7 . 

One can see from this example that accurate estimation 
of third-order moments from turbulence data requires a very 
large record length. Under circumstances where sufficiently 
large data sets are not available, the techniques described 
here and in the last section can be used to estimate the errors 
in the third-moment as quantified by the standard deviation 
03 and the empirical ratio | CT3///3 1 . 

So far in this section the analysis of the third-moment has 
been carried out for one time lag r = 1/50 s. The same 
analysis can be carried out for many different time lags and, 
in this case, the error \cr 3 /fi 3 \ is typically an increasing func- 
tion of time lag t throughout the inertial range (for a fixed 
sample size N). The third-order moment as a function of 
the time lag r computed using all available data is shown in 
the upper plot in Figure 4. For the data shown in Figure 4, 
estimates show that the relative error \a 3 /fi 3 \ lies approxi- 
mately in the range 0.09 < | cx 3 / 1 < 0.3 for r < 0.1 s and 
\a 3 //i 3 \ > 0.3 for r > 0.1 s. Hence, the third-moments are 
sufficiently accurate for the present purpose only for r < 0.1 
s. 

To estimate the energy cascade rate e using Kolmogorov's 
four-fifths law (|2]), the quantity —5((Svh ) 3 )/4r is plotted ver- 
sus t in the lower panel in Figure 4. Note that compress- 
ibility effects are negligible because the Mach number is 
much less than unity and, therefore, the application of Kol- 
mogorov's four-fifths law is justified. In the lower panel in 
Figure 4, the data lie approximately on the horizontal line 
e = 1.8 J/kg-s throughout the range 2 x 10~ 4 < r < 10~ 2 
s. Thus, the value of the energy cascade rate determined by 
Kolmogorov's four-fifths law is e ~ 1.8 J/kg-s. Note that 
the range 2x10~ 4 <t<10~ 2 s where where an apparent 
plateau is formed does not coincide with the inertial range 
10~ 3 < r < 10 _1 s inferred from Figure 1. Because the 
dissipation range lies just beyond the spectral break near 10 3 
Hz in Figure 1 (Pope, page 237), this implies that the region 
where the four-fifths law holds includes part of the dissipa- 
tion range. It is also puzzling why the four-fifths law breaks 
down for r>10~ 2 sin Figure 4 since the inertial range ap- 
pears to extend to r ~ 10 _1 s in Figure 1. Consequently, 
Kolmogorov's four-fifths law does not hold throughout the 
entire inertial range as the theory seems to predict. The rea- 
son for these discrepancies is unknown at the moment. How- 
ever, results for the scaling of the third-order moment in Fig- 
ure 4 are in agreement with lGagne et al.l ( 2004) who analyzed 
the same Modane data. 
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Fig. 2. Histograms of the third-order moments A/3 (AT) for the time lag r = 20 ms obtained from Modane wind tunnel data using samples of 
size N = 2 x 10 3 , N = 2 x 10 4 , N = 2 x 10 5 , and N = 2 x 10 6 . In each case, the total number of moments computed is iV mom , the offset 
between adjacent data blocks is m, and the sample mean and standard deviation /i3 and 03 have units (m/s) 3 . For the case N = 2 x 10 3 the 
range of M3 values in the histogram extends from —70.23 to 49.80 (not shown). The range of M3 values shown in each of the other three 
histograms are all inclusive. 



An independent estimate of the energy cascade rate is ob- 
tained using equation ([T). If the measured signal is approxi- 
mated by the truncated Fourier series 



(N/2)-l 



k=— JV/2 



i2Trkt\ 



(26) 



where N is the record length, T — NAt, and At is the sam- 
pling time, then the time average of (dv/dt) 2 is given by 



q X 2\ /o X 2 N/2 



where 



JV-l 



fe=-(iV/2)+l 



n=0 



N 



(27) 



(28) 



is the discrete Fourier transform of the sequence v n = 
v(nAt) which is easily evaluated using the FFT. Using the 
entire data record to evaluate d27b and the value j/ = 2xlCP 5 
m 2 /s, the energy dissipation rate obtained from equation ([TJ 
is e diss = 2.6 J/kg-s. 

The spatial separation between two consecutive measure- 
ments I = vAt = 0.8 mm is roughly three times the Kol- 
mogorov scale 77 ~ 0.3 mm; the normalized wavenumber 
is krj ~ 2.4. Because most of the dissi pation occurs in the 
wavenumber range krj < 1 dPopell2000l p. 237, Fig. 6.16), 
estimates of ((dv/dt) 2 ) from the Modane data should be ac- 
curate. (Although Fig. 6.16 in Pope's book is drawn for the 
case R\ = 600, a similar plot in the case R\ = 2500 is 
almost indistinguishable from the cas e R\ = 600. ) Note 
that the value of rj given in Table 1 of Malecot et al. ( 200fj|) 
is in error, the correct value is given in Yann Malecot's the- 
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Fig. 3. Empirical results for the quantity 1 0-3/^3 1 as a function of the 
sample size N for the time lag r = 20 ms obtained using Modane 
wind tunnel data (upper plot); the theoretical curve 304/ yN, equa- 
tion (15), is drawn in red. The number AT mom of third-moments 
Ms(N) used to compute the mean /13 and standard deviation 0-3 
are shown in the middle plot. The autocorrelation function of the 
difference series (Svn) 3 is used to determine the correlation time 
t c ~ 90t a used in equation (15); r c is the point where ACF = 0.5 
(bottom plot). 
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Fig. 4. The third-order moment ((<5i>||) 3 } versus time lag r com- 
puted using the entire record of Modane wind tunnel data (upper 
plot). The ratio — 5{(<5«m ) 3 } /4r versus r, where r = vt is the spa- 
tial separation (lower plot). The horizontal line in the lower plot is 
the value 1.8 J/kg-s. The data for r > 0.1 s is uncertain and should 
be disregarded. 



sis and also in IKahalerras et al.l 11998) and in iGagne et al 
(2004) where the same Modane data is used. 



In summary, it has been shown that the energy cascade rate 
e=1.8 J/kg-s obtained by Kolmogorov's four-fifths law is in 
rough agreement with the energy dissipation rate £diss = 2.6 
J/kg-s obtained from equation (Q]). 



6 Illustration using solar wind data 

In this section, we present two examples to illustrate the 
convergence of third-order moments for solar wind data. The 
first example uses 1 second data for the radial magnetic field 
component Br measured above the poles of the sun by the 
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Ulysses spacecraft. The second example uses 64 second data 
for the radial solar wind velocity vr measured in the ecliptic 
plane near 1 AU by the ACE spacecraft. 

6.1 Analysis of Br using Ulysses data 

The radial magnetic field component Br (in spacecraft RTN 
coordinates) was chosen because it enters the third-moment 
(\SB\ 2 SBr) that appears in the law for the cross-helicit 



2007; Podest 



casc ade in MHD turbulence (Podest a et al 
l2008h . Ulysses data was chosen because it is devoid of mag- 
netic sector crossings which are usually present in data ac- 
quired near the ecliptic plane. The third-moment ((5Br) 3 ) 
changes algebraic sign in outward and inward magnetic sec- 
tors. For this reason, the presence of different magnetic sec- 
tors significantly complicates the analysis of this third-order 
moment. 

The Ulysses data selected for analysis consists of a time 
series of ~ 1 s data from the vector-helium magnetometer 



(Balogh et"ail ll992h spanning the time interval from 1 July 
1994 to 1 October 1994, 92 days. During this time Ulysses 
distance from the sun decreased from 2.80 AU to 2.17 AU as 
its heliographic latitude remained between —70 and —81 de- 
grees. The reversal of the solar magn etic field in the southern 
hemisphere was completed in 1992 (ISnodgrass et all 2000) 
so the data used here contains only one magnetic sector. The 
time tags on the data show some data have a 1/2 second ca- 
dence, some data have a 1 or 2 second cadence, and there are 
also data gaps of various sizes. The 1/2 sec data is downsam- 
pled to 1 sec and the data gaps are left intact to create a time 
series with a uniform cadence of 1 sec. Times when data 
are missing are marked with fill values (such data are easily 
omitted from the analysis). There are a total of 7.95 x 10 6 
data points in the time series and 23.55% of these points are 
missing data markers (fill values). The average value of Br 
for the entire time series is —0.526 nT. 

The power spectrum for the Ulysses data shown in Figure 
5 is strikingly similar to that of the wind tunnel data in Figure 
1. From Figure 5, the inertial range appears to extend from 
less than 10~ 3 Hz to approximately 10 _1 Hz. The time lag 
t — 60 seconds or, equivalently, / = r _1 = 1.67 x 10 -2 Hz 
is chosen for analysis because it lies inside the inertial range. 
The same procedures used to analyze the Modane data are 
employed for the Ulysses data except that missing data is ex- 
cluded from the analysis. Consequently, a data block of size 
M contains less than M samples (because of the presence of 
fill values) and the actual number of samples N varies from 
block to block. Only those data blocks where N > 0.55M 
are included in the analysis and the average number of sam- 
ples N taken over all blocks of a given size M is defined to 
be the sample size N for that run. 

The results of the statistical analysis of Ulysses data for the 
time lag t = 60 seconds are shown in Figure 6. The sizes 
of the data blocks used in the analysis are M = 2 x 10 4 , 
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Fig. 5. Power spectrum of the radial magnetic field component Br 
for the Ulysses data used in this study (blue dots). The red line is 
proportional to /~ 5 ^ 3 . The inertial range extends from less than 
10~ 3 Hz to approximately 10 _1 Hz. The onset of the dissipation 
range is indicated by the change in slope around 10" 1 Hz. 



2 x 10 5 , and 2 x 10 6 . The offset from one data block to 
the next is M for the upper plot and Mj 2 for the other two 
plots. The algebraic sign of /i3 is negative except in the lower 
plot, however, the sample size in the lower plot is too small 
to yield adequate statistics. The value of the third-moment 
obtained using the entire data record is —8.9 x 10~ 5 (nT) 3 , 
a very small value. To gain some idea of the error, the error 

1 /2 

of the mean (Ta/Nmom estimated from Figure 6 is roughly 
2 x 10~ 4 (nT) 3 . 



The theoretical relation (15) may be evaluated using es- 
timates obtained from the data for the third-moment M3 = 
-8.9 x 10~ 5 (nT) 3 , the sixth-moment M 6 = 2.1 x 10~ 2 
(nT) 6 , and the correlation time r c ~ 21s determined from 
the autocorrelation function for the sequence (5 Br) 3 . The 
values of the moments are uncertain, especially higher order 
moments such as M 6 which ca n be strongly affected by th e 
presence of outliers in the data dHorburv and Balo gh. 1997), 
however, they are used anyway to explore the fit to the data 
of the relation (15). Thus, the theoretical relation (15) takes 
the form 



<7 3 (AO 



MAO 



7500 



N 



(29) 



For the three runs shown in Figure 6 the theoretical values 
for the ratio |cr 3 / // 3 1 are 57, 19, and 6. These are roughly 
consistent with the values found in Figure 6. For all data, 
N ~ 6 x 10 6 and the relation d29l predicts a relative errror 



J. J. Podesta: Accurate estimation of third-order moments from turbulence measurements 



9 



70 
60 
50 



g 40 
8 30 



20 
10 




15 



N = 1.73x10' 
|a 3 /u 3 |=25 



'3 

mom 



5.3x10 
282 



-0.02 




0.02 



M 3 (nT) 



10 



c 

o 
O 



N = 1 .52x10° 

l*Al = 17 4 

H 3 =-1-2x1 0" 4 
a =2.1x1 0" 3 

NT =78 

mom 



I 



-10 



M 3 (nT) J 



x 10" 



1.5 



c 

=3 

o 
O 



0.5 



N = 1.49x10° 
|a 3 /u 3 | = 2.1 

H =2.4x1 0" 4 

a, = 5.0x1 0" 4 

Nr =6 

mom 



-0.5 



0.5 



M 3 (nT) J 



x 10" 



Fig. 6. Results from the analysis of the third order moment 
((5Br) s ) for the time lag r = 60 s obtained using data from 
Ulysses first southern polar pass. The number TV is the approxi- 
mate sample size used to compute the third-moment and iV mom is 
the number of third-order moments used to compute the statistics 
pi3 and 03. The mean and standard deviation of A4s(N), denoted 
by Hz and 03, have units (nT) 3 . The value of M3 obtained using the 
entire data record is -8.9 x 10~ 5 (nT) 3 . 



\a 3 /fi 3 \ ~ 3. Hence, fi 3 ^ -8.9 x 1CT 5 ±2.7 x 1CT 4 (nT) 3 , 
which is consistent with the estimates in the preceeding para- 
graph. 



In summary, the large magnitude of the ratio (J3//X3 is 
partly due to the fact that /13 is close to zero and this makes it 
impossible to obtain adequate convergence with the limited 
data used in this study. One may conclude from these results 
that a much larger data set than the one used here is needed 
to determine the third-moment of 5Br accurately. Neverthe- 
less, the results presented here are still useful for determining 
approximate upper and lower bounds for this third-order mo- 
ment. 

6.2 Analysis of vr using ACE data 

Solar wind measurements of the radial velocity component 
from the Advanced Composition Explorer (ACE) are ana- 
lyzed in the same way. The ACE spacecraft is in orbit around 
the Sun-Earth L\ libration point 240i? e sunward of the Earth. 
The A CE SWEPAM plasma i nstrument has a 64-second ca- 
dence dMcComas et all [l998) and we use all data available 
during the three year period from 2005 through 2007, about 
1.4 million data points. Note that solar minimum is expected 
to occur in late 2008 or early 2009. Non-overlapping data 
blocks of 100 points (about 10,000 blocks) to 256,000 data 
points (5 blocks) are used. Each data block may include fill 
values (missing data markers) that are present in the time se- 
ries. All fill values are omitted from the analysis and any 
data block in which the number of fill values exceeds 10% 
is excluded from analysis. Third-order moments of 5vr are 
calculated for two different time lags, r = 256 seconds and 
r = 2048 seconds. The inertial range in the ecliptic plane 
near 1 AU extends from about 1 second to about 1 hour and, 
therefore, both of these time lags lie in the inertial range. 

Figure 7 shows that in the solar wind, the ratio of the stan- 
dard deviation of the third-moment to the average value of 
the third-moment \a 3 /fj, 3 \ has the same N^ 1 / 2 dependence 
predicted by equation (15) as does the wind tunnel data an- 
alyzed in Section 5. Remarkably, the amplitude of this re- 
lation is quantitatively similar for both wind tunnel data and 
solar wind data, even though the solar wind has a fast/slow 
stream structure and the turbulence is magnetohydrodynamic 
in nature. As with the wind tunnel data, it appears from 
Figure 7 that more than 10 7 solar wind velocity measure- 
ments are needed to accurately determine the third-order mo- 
ment ({Svr) 3 ) (errorless than 10%). However, sample sizes 
N ~ 10 6 may give sufficient accuracy for some applications. 



Also plotted in the upper plot in figure 7 are the theoret- 
ical curves, equation (15), for the two time lags studied. A 
rough estimate of equation (15) is obtained by using all data 
in the record to estimate the sixth-moment of Svr(t), Mq, 
the third-moment M 3 , and the correlation time t c of the se- 
quence [Svh(t)] 3 . For the time lag t — 256 s, this yields 
M e ~ 8.8 x 10 9 (km/s) 6 , M 3 ~ 2.5 x 10 2 (km/s) 3 , and 
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Fig. 7. Results from the analysis of the third-order moment 
({Svr) 3 } using ACE data for 2005-2007. Red triangles and blue 
squares correspond to the time lags r = 256 s and r = 2048 
s, respectively. The upper plot shows the convergence of the ratio 
|cr3 //Lt3 1 as a function of the number of samples N; the solid lines 
are the theoretical predictions from equation (15), 477/A 1//2 and 
322/ TV 1 / 2 . The lower plot shows the values of the third-order mo- 
ment plus and minus the error in the mean, that is, fis ± 0-3 /JV™, 
where A mom is the number of moments used to compute the mean. 



t c /t ~ 0.4. In this case, n — 1.6 and equation (15) becomes 



M3(A0 



477 



(30) 



For the time lag r = 2048 s, M 6 ~ 9.5 x 10 10 (km/s) 6 , 
M 3 ~ 3.4 x 10 3 (km/s) 3 , and t c /t ~ 0.4. In this case, 



12.8 and equation (15) becomes 



322 

Vn' 



(31) 



These represent reasonable asymptotic fits to the data shown 
in Figure 7 as A becomes large. 

What accuracy can be claimed for measurements of this 
third-order moment in the solar wind? The upper plot in Fig- 
ure 7 indicates that for N ~ 1.4 x 10 6 data points the error is 



around ~ 40%. The lower plot in Figure 7 shows the mean 
values of the third moment /13 plus and minus the standard 

1 /2 

error of the mean (73/AV om , where N mom is the number of 
moments used to compute the mean. It appears that the rel- 
ative error is roughly the same at both lags and that within 
the calculated errors the third moment is proportional to lag. 
Such proportionality is most clearly demonstrated by com- 
puting the third-moment as a function of the time lag r (not 
shown). 

It is interesting that both R\ ~ 2500 wind tunnel data and 
solar wind velocity data require roughly the same number of 
data points to obtain good convergence of third-order mo- 
ments for time lags in the inertial range. In part, this may be 
because both kinds of turbulence have similar Reynolds num- 
bers. The Reynolds number in the solar wind can crudely be 
estimated using the hydrodynamic relation L/r/ = Re 3 / 4 , 
where L is the integral scale, 77 is the Kolmogorv scale (dis- 
sipation scale), and Re is the Reynolds number based on 
the integral scale. Solar wind power spectra indicate that 
L/rj ~ 10 5 and, therefore the Reynolds number is of order 
10 6 . This is equivalent to a Taylor-scale Reynolds number 



Rx = (fRe 



.1/2 



2600 



(32) 



(Pope, 2000l p. 200, eqn. 6.64). Thus, the Reynolds numbers 
of Modane wind tunnel data and solar wind velocity data at 
1 AU are similar. 



7 Conclusions 

The purpose of this study is not to compute turbulent en- 
ergy cascade rates using third-order moments. The purpose 
of this study is to show how the accuracy of third-order mo- 
ments can be estimated and controlled. A simple theory is 
presented that describes the statistical convergence of third- 
order moments, such as ((Svn ) 3 ), as a function of the record 
length. An important conclusion is that the accuracy of third- 
order moments depends on the number of correlation lengths 
spanned by the time series as expressed by equations (14) 
and (15). The techniques described here are useful for as- 
sessing the accuracy of third-order moments obtained using 
measured data. Practical applications of the theory have been 
illustrated using wind tunnel data and solar wind data. 

For the accurate calculation of third-order moments from 
wind tunnel data with a Taylor-scale Reynolds number R\ ~ 
2500, the number of data points required to obtain an error 
less than 10% at the time lag r = 20 ms is N = 10 7 or, 
equivalently, a record length spanning N c = 10 5 correlation 
lengths. For fluctuations of the radial solar wind velocity 
Vr, the analysis of ACE data in the ecliptic plane near 1 AU 
shows that for the time lags r = 256 s and r = 2048 s the 
number of data points required for an accurate determination 
of the third-order moment is also roughly A" = 10 7 . This 
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is equivalent to N c = 6.3 x 10 6 and 7.8 x 10 5 correlation 
lengths for r = 256 s and 2048 s, respectively. However, 
N ~ 10 6 data points may yield sufficient accuracy for some 
applications. 

For fluctuations of the radial magnetic field component 
Br over the poles of the sun at a heliocentric distance of 
approximately 2.5 AU, the value of the third-order moment 
is close enough to zero that convergence of the third moment 
could not be demonstrated using an interval of Ulysses data 
with approximately six million points (not including fill val- 
ues), a record consisting of approximately 3 x 10 5 correla- 
tion lengths. This suggests that third-order moments of solar 
wind magnetic field components must be computed carefully 
because without a sufficiently large number of data points 
and without evaluation of the probable errors using eq (15) 
the calculation of the third-order moments are not meaning- 
ful. 

It should be noted that the above two examples based on 
solar wind data from Ulysses and ACE are distinctly differ- 
ent from each other and from the example based on Modane 
wind tunnel data. The Ulysses study pertains to the radial 
magnetic field component at high heliographic latitudes and 
the ACE study pertains to the radial velocity component in 
the ecliptic plane. These studies provide two separate ex- 
amples of the estimation of third-order moments and their 
uncertainties using solar wind data. For wind tunnel data, 
Kolmogorov's four-fifths law predicts the third-order veloc- 
ity moment scales linearly in the inertial range. For solar 
wind data, neither the third-order moment of the radial ve- 
locity component nor the third-order moment of the radial 
magnetic field component is predicted to scale linearly. The 
latter examples simply serve to illustrate the application of 
statistical convergence techniques to solar wind data. It is 
of interest to note, however, that the third-order moment of 
the radial velocity component in the solar wind was found to 
scale approximately linearly in the study by MacBride et al. 
(2008) and appears to provide the dominant contribution to 
the energy cascade rate estimated from the scaling relations 
of Politano and Pouquet (1998). 

In conclusion, the present study has some noteworthy im- 
plications for measurements of the energy cascade rate in the 
solar wind. Empirical estimates of the energy cascade rate in 
the solar wind have recently been obtained under the assump- 
tion that the turbulence is approximately incompressible and 
isotropic so that the third-moment scaling relations of Poli- 
tano and Pouquet (1998) for homogeneous isotropic incom- 
pressible MHD turbulence cou ld be applied. In the studies by 
Sorriso-Valvo et all d2007b and lMarino et al.1 d2008l) the num- 
ber of data points employed to compute the required third- 
order moments was around 2000. As shown in the present 
study, this number is insufficient to obtain accurate estimates 
of thi rd-order moments in the s olar wind. This may e xplain 
whv lSorriso-Valvo et all d2007t) and lMarino et alj d2008l) did 



not find a linear scaling of the third-order moments in some 
of the intervals they studied, and why they found it was rare 
for linear scaling to be observed simultaneously for both of 
the Elsasser variables. Although MacBride et al. (2008) did 
not use the convergence tests proposed here, they used large 
enough data sets that the third order moments in the Poli- 
tano and Pouquet scaling laws became insensitive to adding 
more data. To obtain stable estimates of the third-moments, 
this convergence criterion required the use of at least one 
year of ACE plasma and magnetic field data, roughly 10 6 
data points. In the future, the convergence of third-order mo- 
ments and the associated error estimates that such conver- 
gence studies provide should become an integral part of any 
analysis of solar wind data involving third-order moments. 
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